Energy diffusion in strongly driven quantum chaotic systems: 
Role of correlations of the matrix elements 

P.V. Elyutin and A.N. RubtscJ] 

Department of Physics, Moscow State University, Moscow 119992, Russia 

Vq- (Dated: February 8, 2008) 

O ' 

O ■ Abstract 

O ■ The energy evolution of a quantum chaotic system under the perturbation that harmonically 

<D ■ 

Q , depends on time is studied for the case of large perturbation, in which the rate of transition 

00 ' 

calculated from the Fermi golden rule (FGR) is about or exceeds the frequency of perturbation. For 
this case the models of Hamiltonian with random non-correlated matrix elements demonstrate that 
the energy evolution retains its diffusive character, but the rate of diffusion increases slower than 
the square of the magnitude of perturbation, thus destroying the quantum-classical correspondence 
for the energy diffusion and the energy absorption in the classical limit H — * 0. The numerical 
calculation carried out for a model built from the first principles (the quantum analog of the 
fSJ , Pullen - Edmonds oscillator) demonstrates that the evolving energy distribution, apart from the 

diffusive component, contains a ballistic one with the energy dispersion that is proportional to 

(N ' 

the square of time. This component originates from the chains of matrix elements with correlated 
signs and vanishes if the signs of matrix elements are randomized. The presence of the ballistic 
component formally extends the applicability of the FGR to the non-perturbative domain and 
restores the quantum-classical correspondence. 
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I. INTRODUCTION 

The problem of susceptibility of chaotic 
systems to perturbations has attracted much 
attention in the last decade [1 - 9]. This 
problem is fundamental, since it includes the 
determination of the response of a mate- 
rial system to an imposed external electro- 



magnetic field, the setup that is typical for 
many experiments. Due to the sensitivity of 
classical phase trajectories or quantum en- 
ergy spectra and stationary wavefunctions of 
chaotic systems to small changes of their pa- 
rameters, the problem is challengingly diffi- 
cult. A consistent and noncontroversial pic- 
ture covering (albeit qualitatively) all the es- 
sential cases of the response hasn't been yet 
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drawn at present. From the point of view ergy increases linearly with time, {AE 2 (t)) = 

of general theory, the problem is related to 2Dt, where D (E, F, uj) is the energy diffusion 

the applicability of the concept of quantum- coefficient. 

classical correspondence to chaotic systems, If the external field F is sufficiently small 

that is a long-standing question in its own in comparison with the appropriately aver- 



right [lCj, LL3J. aged values of the forces acting on a particle 



We shall study a one-particle system with in the unperturbed system, then in the clas- 

the Hamiltonian of the form H = H - sical model the eneT ^ diffusion coefficient D 

Fxcosut, where F (p,f) is the Hamilto- can be expressed through the characteristics 

nian of the unperturbed system; p and f are of the unperturbed chaotic motion of the ac- 

the operators of Cartesian components of the ^ ve coordinate, namely 

momentum and of the position of the par- n „ , 

D = -lu 2 F 2 S x (E,u), (1) 

ticle. The classical system with the Hamil- 
tonian function Hq (p, r) will be assumed to where S x (E,u) is the power spectrum of 
be strongly chaotic, that is, nearly ergodic the active coordinate (the Fourier transform 
on the energy surfaces in a wide range of of its autocorrelation function) for the mo- 
the energy values, system with d > 2 de- tion on the surface with the constant energy 
grees of freedom. In the perturbation opera- value E [9]. The same expression (1) in the 
tor V (t) = —Fx cos out the active variable x case of weak perturbation can be obtained in 
is one of the Cartesian coordinates of the par- the classical limit from the quantum model, 
ticle, coupled to the external homogeneous The evolution of the quantum system can be 
force field. The amplitude F in the following treated as a sequence of one-photon transi- 
will be referred to as the field. In the fol- tions between stationary states of the unper- 
lowing we shall deal with the quasiclassical turbed system |n) — > \k), accompanied with 
case, when the Planck constant h is small in the absorption or emission of the quanta hu. 
comparison of the action scale of the system For small h the energy spectrum of Hq is qua- 
Hq. sicontinuous, thence the rates of transition 

Under the influence of the perturbation are S iven b ^ the Fermi S olden mle ( FGR ) 

the energy value E(t) = H (t) varies in a • tt 2 2 x /n x 

W F = tttF \x nk \ p{E k ), (2) 

quasirandom way. These variations can be 

frequently described as a process of the en- where x nk is the matrix element of the ac- 

ergy diffusion |l2,[l^|. when for the ensemble tive coordinate, and p{E k ) is the density 

with the microcanonical initial energy distri- of states near the final state of the tran- 

bution H (0) = E the dispersion of the en- sition. Although the matrix elements x nk 



in quantum chaotic systems fluctuate wildly Weisskopf - Wigner formula A = hW 16|, 



with the variation of k [10, 11], the averaged is small in comparison with the quanta en- 

squared quantity \x nk \ in the limit h — > is ergy hu. From Eq. (4) it is evident that in 

smooth; it is proportional to the power spec- the classical limit h — > this condition will 

nn 

trum S x (E, u) of the coordinate [fj, H5| , be violated. In the following we shall use the 

S x (E,u) border value of the field F b , defined by the 



ftp (-£/) condition (Ft,) = uj , and refer to the do- 

From Eqs. (2) and (3) we have for the tran- main F > F b as the range of the strong field. 

sition rate By analogy with other models, for strong 

71- „ fields one can expect a slow-down of the 
W F = —^F 2 S x (E,u). (4) 

2n growth of the energy diffusion coefficient D 

Then for the energy dispersion for small t we and of the energy absorption rate Q. For 

have (AE 2 ) = 2D F t = 2 (hu) 2 W F t, that re- example, for a two-level system with relax- 

turns us to the Eq. (1) for the energy dif- ation the quadratic dependence the absorp- 

fusion coefficient. It can be shown that the tion rate Q oc F 2 for small field turns into 

same expression for D holds also for large t a field-independent value Q for strong one. 

The energy absorption in chaotic systems The border is determined by the condition 

comes as an epiphenomenon of the energy dif- £l 2 /TiT 2 ~ 1, where Q is the Rabi frequency 

fusion [4]. With the account of the depen- and Fi,r 2 are longitudinal and transversal 

dence on the energy of the power spectrum relaxation rates correspondingly The 

S x (E, u) and the density of states p (E) the rate of transitions from discrete to continu- 

diffusion becomes biased, and the energy ab- ous energy spectrum (that are basically co- 

sorption rate Q is given by the formula jg], [4] variant with the energy absorption rate Q), 

1 j studied in the context of the theory of pho- 

^ p dE ^ ^ ^ toionization, for sufficiently strong fields can 

Although for weak fields D does not de- even decrease with the increase of F - the ef- 

pend on the Planck constant h, the condi- fect that is known as atom stabilization by 

tion of the applicability of Eq. (2) does, the strong field [lj. 

The FGR is, after all, only a formula of the The slow-down of the energy diffusion in 

first order perturbation theory. It is based strong harmonic fields for the model of quan- 

on the assumption that the transition pro- turn chaotic systems with random uncorre- 

cess has a resonant character - that the width lated matrix elements has been first demon- 

A of the energy distribution of states pop- strated by Cohen and Kottos [5|. A different 



ulated from the original one, given by the approach 



19] has lead to qualitatively the 



same results. This slow-down destroys the sure of the chaotic component ^ S {E) on the 

quantum-classical correspondence. surface of Poincare section) and in intensive 

It has been demonstrated by Kottos and (that is measured by the magnitude of the 

Cohen 20| that for the first-principles model Lyapunov exponent cr{E)) aspects. For val- 

that is constructed by the quantization of the ues of energy E > 2.1 the measure fj, s > 0.5, 

Hamiltonian of a classically chaotic system, and chaos dominates in the phase space; for 

the response to a sudden change of the (oth- E > 5 the chaotic motion of the system is 



22J. 



erwise) stationary Hamiltonian measured by approximately ergodic 
the energy spreading restores its classical be- The matrix of the quantum Hamiltonian 
haviour for sufficiently small values of H in operator of the model Eq. (6) has been cal- 
contrast with the model with random inde- culated in the basis of the unperturbed two- 
pendent matrix elements. dimensional isotropic harmonic oscillator for 

The purpose of the present article is to the value H = 0.05. Due to the symmetry 

study this phenomenon for the harmonically of the system the submatrices with different 

driven system. parities of the quantum numbers n x and n y 

can be diagonalized separately. 

II. NUMERICAL EXPERIMENT By expanding the wavefunction of the sys- 
tem \1/ (r, t) in the basis of the eigenstates 

The system chosen is the Pullen - Ed- {^j of the system # q 

monds oscillator 21(, that describes the two- 



dimensional motion of a particle in the quar- * ( r ' X! a m(t)(pm (r) e m , (7) 

m 

tic potential. The Hamiltonian of this system 



is 



we obtain for the amplitudes a m the system 
of equations 

H ° = L {&+fi) + ^T { x2 + y2 + X ~W) ■ ^^ = £lWoso;^a m , (8) 

(6) k 

In the following we use the particle mass m, where the quantities flkn — ^ X Fx^ n are the 
the frequency of small oscillations cuo, and the Rabi frequencies of transitions. This system 
nonlinearity length A as unit scales, and write has been solved numerically for the 2352 am- 
all equations in dimensionless form. plitudes of the eigenstates with "even-even" 

The properties of chaotic motion of the and "odd-even" parities of n x and n y that in- 
Pullen - Edmonds model are thoroughly stud- elude all states of these classes with energies 



ied 



22 



23, 



24l | . With the increase of energy in the band 10 < E < 12. 



the system becomes more chaotic both in ex- For the initial conditions in the runs with 
tensive (that is characterized by the mea- different values of F and u we have used 



the same normalized narrow wavepacket with 
randomly chosen real amplitudes Ofc(0) with 
the complexity (inverse participation ratio) 
C = (j^aj. (0)J = 21, the mean initial en- 
ergy (E) = E = 11.0 and the initial energy 
dispersion AE 2 (0) = 2.5 x 1(T 3 . 

To expose the role of correlations of the 
matrix elements Xkn on the energy kinetics 
we compare the properties of the "natural" 
system with its "randomized" analog with 
the matrix elements y mn = x mn A mn , where 
A mn are the elements of a symmetric matrix 
that take values 1 or -1 at random with equal 
probabilities. This randomization, that has 
been introduced in {20], destroys the correla- 
tion of matrix elements. The tests has shown 
that the results do not depend on the specific 
choice of the matrix A mn within the limits of 
an error of about 1%. 

Figure 1 presents the typical dependence 
on time of the energy dispersion for the nat- 
ural (filled circles) and randomized (open cir- 
cles) models. To extract the value of the en- 
ergy diffusion coefficient, the numerical de- 
pendence was fitted by the law of evolution of 
the energy dispersion for the diffusion equa- 
tion with the constant D on the interval 
(— L, L) with impenetrable walls on the bor- 
ders and the initial condition in the form of 
the 5(0) peak. This dependence can be de- 
scribed (with the local accuracy better than 
3%) by the formula 

$(a,b;t) = a(l -exp(-fot)) x (9) 
(1 + 0.6336* exp (-1.1616*)) 
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FIG. 1: The dependence of the energy disper- 
sion AE 2 on time t for the natural (filled cir- 
cles) and randomized (open circles) models. The 
energy dispersion is measured in the units of 
(Tluiq) 2 . The field strength corresponds to the 
Fermi point: F = Fj, = 0.022; the perturbation 
frequency uj = 1.00. 

where a = L 2 /3 and 6 = 6D/L 2 . The best 
fits of the formula Eq. (9) with the numerical 
data are shown in the Fig. 1 by solid lines. 
The initial moment *o has been used as the 
third fitting parameter. 

For once, it is clearly seen that the ran- 
domization suppresses the process of the en- 
ergy diffusion. Secondly, for the natural 
model one can see the presence of two differ- 
ent regimes - a fast initial diffusion sharply 
slows down at a crossover time t c ~ 20. We 
note that at this moment the energy disper- 
sion AE 2 (t c ) = 32 (TiuJo) 2 is much less than 
the saturated value AE 2 = a = 133 (Hlvo) 2 , 
that corresponds to the uniform probabil- 
ity distribution throughout the band of the 
states taken into account. To understand 



what happens at the crossover time we have 
to study more closely the time development 
of the probability distribution. 

The overall form of the energy distribu- 
tion for the natural model is rather accurately 
approximated by the gaussian form that fol- 
lows from the model of the energy diffusion 
with the constant D. This agreement can be 
seen in Fig. 2, where the logarithm of the 
energy density distribution is shown against 
the reduced energy shift Ae — (E — E ) /hcuo, 
where E is the mean energy value of the ini- 
tial wavepacket. 

* °l 




Reduced energy shift Ae 

FIG. 2: The logarithm of the energy density dis- 
tribution for the time t = 17.3 for the parameters 
F = 0.606F&, ijj = 1.00. The grassy line is the nu- 
merical data, the solid line presents the best fit 
of data with the parabola p(Ae) = a — (3{Ae) 2 . 

Although the agreement seems to be very 
good, one must keep in mind that the ver- 
tical scale of the graph is logarithmic. By 
subtraction of the parabolic fit from the nu- 
merical distribution we come to the picture 
of deviations that is shown in the Fig. 3. 



■ nil 


1 ,JI 


i, lli, 




If! 1 "! 




If 



-20 ^10 10 20 

Reduced energy shift Ac 



FIG. 3: The deviation of the logarithm of the 
energy density distribution from the best fitted 
parabola Aw as a function of the dimensionless 
energy deviation Ae for the same parameters as 
in Fig. 2. The black stars note the maxima of the 
ballistic bumps (see the text) 

The peak at Ae = corresponds to the 
part of the initial packet that is not depleted 
yet by the energy spreading. Two bumps are 
clearly seen in the picture: their maxima are 
located at Ae ~ ±12. The calculations show 
that these bumps propagate outwards with 
a constant speed, that is ballistically. The 
crossover time corresponds to the moment 
when these bumps reach the borders of the 
treated band of states. Therefore, only the 
part of the graph that precedes the t c cor- 
responds to the properties of the model; the 
second part is just an artifact of the calcula- 
tion scheme. To estimate the diffusion coef- 
ficient we fitted the dependence of the initial 
state by the two-parameter formula 

2Dt 2 

4>(D,T-t) = — (10) 

T + t 



where the time shift r accounts for the dura- equidistant energy spectrum E n = Tiujq [n/2] 

tion of the initial stage, where the law of the where [, ] denotes the integer part of the num- 

dispersion growth is always quadratic. The ber. The matrix elements of the coordinate 

best fitted function is plotted in Fig. 1 by connect each state \n > to all four states 

the dashed line. \m > with the energy differences E n — E m = 

The ballistic spreading of the energy ±hu : 
distribution in the time domain t » fx x x x \ 

lixjuj is well known for the model of a , , 

\\6) 

one-dimensional resonantly excited harmonic „ , 

tor even n and 

oscillator with the Hamiltonian H (t) = 

(f + x 2 ) /2 - Fx cos t . In the quasiclas- Xnm = X ^ n ' n ~ 3 + 5n > n ~ 2 + 5n < n + l + <W+ 2 ) 

sical domain (for large quantum numbers 

,\ ,, , • , , r , i for odd n. In this model for the resonant 

n 3> 1) the matrix elements of the co- 

,. , i i i , , perturbation V (t) = —Fx cos uint the energy 

ordinate can be taken constant, x nm — ^ w u toJ 

v , r . r \ , c ' 1 1 , spreading is ballistic, AE 2 (t) oc t 2 , whereas 

X (o n n -i + On,n+i)> where dij is the Kroneker K to ' w ' 

, i , ttt-,1 ,i • . . the randomization of signs leads to the lo- 

delta - symbol. With this assumption m the 

, . . . r , i . . , . , calization of the quasienergy states, and the 

rotating wave approximation tor the initial 

j., /n x r i . energy dispersion growth saturates, 

condition (U) = Ofc n the probabilities Wk to OJ 1 ° 

c j , i , ii ii • u O ne can suppose that the ballistic compo- 

find the system m the state \k > are given by 11 v 

, , i, i r , nent in the quantum chaotic model is carried 

the well-known formula n 

through the subset of states that are similar 
w k = \a k {t)\ 2 = Jl_ k {Slt) , (11) to the "double ladder" model. The degree 

of correlation of the matrix elements can be 

where J n (z) are the Bessel functions of the 
first kind and f2 = FX/2K is the value of the 
Rabi frequency. Equation (11) yields for the = %nm%miXikXkn, (15) 

m,l,k 



estimated from the construction 



energy dispersion 



that describes the sum of contributions of 



^£>2 ^ _ 1 {fuj^t) 2 (12) a ^ P oss ible four consequent transitions that 

start and end on the same state \n >. The 
However, the randomization of signs of ma- corre i at i on in dex v can be defined as the ratio 
trix elements does not influence the energy of the average value of T for the ran domized 
kinetics in this model. system to that of the natural system. For the 

The influence of randomization can be ex- "double ladder" model we have 

plained by the toy "double ladder" model. (t r \ 7 

v = n = — = 583 (16) 
This system has the doubly degenerate (^ v ) 12 



The relatively large value of this number is 
explained by the large contribution of sym- 
metric contours with m = k, that are invari- 
ant under the randomization. For the Pullen 
- Edmonds model the value of the correlation 
index v = 0.72 is rather close to that of the 
"double ladder" model; that makes the anal- 
ogy plausible. 

It must be stressed that the combined ef- 
fect of the bulk diffusion spreading and the 
overlaying ballistic bumps propagation pro- 
duces the linear growth of the energy disper- 
sion (see Fig.l) that will be referred to as the 
effective diffusion. 

We define the repression coefficient 
R(F, uj) as the ratio of the effective energy 
diffusion coefficient to its value Dp that fol- 
lows from the FGR: 
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FIG. 4: The dependence of the repression coef- 
ficient R on the logarithm of the field strength 
A = ]n(F/Fb) for the frequency of perturbation 
u = 1.00. 



R (F, uo) 



D(F,u) 



(17) 



2(huyw F 

The dependence of the R on the field strength 
is shown in Figs. 4 and 5 for two different 
values of the perturbation frequency. 

For the natural model R remains approxi- 
mately constant with the value close to unity, 
whereas for the randomized model the energy 
diffusion slows down in qualitative agreement 
with the conclusions of p, Q . 

III. ENERGY DIFFUSION IN THE 
CLASSICAL MODEL 

The classical expression for the diffusion 
coefficient Eq. (1) is derived in the limit of 
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FIG. 5: The dependence of the repression coef- 
ficient R on the logarithm of the field strength 
A = ]n(F/Fb) for the frequency of perturbation 
u = 1.62. 

the infinitesimal perturbation, when one can 
neglect the influence of the perturbation on 
the law of motion of the active coordinate 
x(t). Let's study the formation of this coeffi- 
cient. We represent the external field in the 
form Fsinut and denote by t n = 2-nnjuj the 
moments of time at which the external field 



take zero values. The variation of the energy 
for one field period between these moments 
is exactly proportional to the field strength, 

t n +l 

AE n = F J x (t) sin utdt = FA n . (18) 

The quantities A n we shall call the reduced 
variations of the energy. In the accepted ap- 
proximation they do not depend on the field 
strength. 

The following calculations were carried 
out for the Pullen - Edmonds oscillator on 
the energy surface E = 11.0 and for the per- 
turbation frequency u = 1.00. 
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FIG. 6: The dependence of the autocorrela- 
tor of the reduced energy increments (k) = 
(A n A n+ fc) on the time shift k measured in pe- 
riods of the perturbation on the energy sur- 
face E = 11.0 for the perturbation frequency 
co = 1.00. 

The values of the reduced energy incre- 
ments A n on the neighbouring time intervals 
are correlated. Figure 6 presents the form of 
the autocorrelation function of the reduced 
energy increments. One can see that for the 



values of the time shift k > 5 the correlations 
become rather small. 
The quantity 

<* = ^(t*~<)' < 19 > 

will be called the K-th approximant of the 
reduced diffusion coefficient. This is a pro- 
portionality coefficient between the diffusion 
coefficient and the square of the field ampli- 
tude, calculated from the interval of time of 
K consequent periods of the field. The pos- 
itive correlation of A& for small k produces 
the initial monotonous growth of the that 
rather rapidly comes to a saturation. 
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FIG. 7: The dependence of the approximants of 
the reduced diffusion coefficient d n on the du- 
ration of the time interval measured in periods 
of the perturbation for the Pullen - Edmonds 
model Eq. (6) on the energy surface E = 11.0 
for the perturbation frequency uj = 1.00. 

From the graph in Fig. (7) it is seen that 
already ds takes the value that within the 
1.5% error margin is undistinguishable from 
the asymptotic limit. However, this average 



quantity is formed by the contributions that 
differ by several orders of magnitude. The 
graph in Fig. 8 shows the distribution of the 
quantities d§ in the log-log scale. The dis- 
tribution is taken from averaging over four 
ensembles of 10 5 points each. 
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FIG. 8: The distribution of the values of partial 
contributions to the approximant ds from differ- 
ent points on the trajectory of the Pullen - Ed- 
monds model on the energy surface E = 11.0 for 
the perturbation frequency u = 1.00. The thin 
lines show the approximate forms of the distri- 
bution given by the formulas in the text. The 
black star marks the maximal partial contribu- 
tion <i + that comes from the resonant trajectory, 
the white star marks the position of the average 
value < ds >■ 

The dominating part of this distribution 
is accurately fitted by the dependence In w — 
/(d) = -l.771-l.0571nd-0.095 (hid) 2 that 
is shown by the thin solid line. For the largest 
values of d$ another approximation is valid, 
\nw = g(d) = 19.268 - 6.386 In d. This 
dependence is shown in Fig.8 with the thin 



dashed line. In the domain of validity of the 
approximation /(d) the slope of the curve is 
less than unity: the distribution is of the Zipf 
- Pareto type, in which the dominating con- 
tribution to the average comes from the rare 
large terms. In our case the 20% of the largest 
terms come with 82% contribution to the av- 
erage. These large contributions come from 
the bits of the trajectories in which the point 
oscillates almost along the direction of the 
perturbing force nearly synchronously with 
the perturbation. Theoretically the maxi- 
mal value of ds originates from the motion 
with the law x(t) = \f2E sin(t) and equals to 
d+ = AixE / uj = 138.23. 

In this resonant case the energy increment 
grows linearly in time - that is, ballistically. 

Thus we can indicate a classical coun- 
terpart to the quantum dynamics of energy 
growth. The quantum ballistic bumps are 
analogous to the nearly resonant bits of the 
classical trajectories with the quasiballistic 
energy increase. 

IV. CONCLUSION 

By the numerical studies of the evolution 
of the energy distribution in a harmonic ex- 
ternal field in a system constructed by the 
quantization of a classically chaotic Hamil- 
tonian system, thus retaining all correlations 
of the matrix elements, we have found that 
the effective rate of the energy diffusion pre- 
serve its quadratical dependence on the field 



strength on the domain of the strong field, 
where the transition rate is comparable to the 
perturbation frequency. In other words, the 
Fermi golden rule appears to be valid far be- 
yond the limits of the domain in which its 
applicability can be justified. This circum- 
stance restores the quantum - classical cor- 
respondence for the energy diffusion and the 
energy absorption rate in the limit h — > 0. 

We have to admit that our studies are lim- 
ited to a specific model, studied for two val- 
ues of the perturbation frequency. However 
the revealed mechanism of the ballistic com- 
ponent of the energy distribution that prop- 



agates through the chains of matrix elements 
with the correlated signs can be admitted as a 
universal, especially with the account of the 
important contribution of the quasiballistic 
parts of the trajectories to the energy diffu- 
sion in the classical model. The more rigor- 
ous proof of the universality demands further 
studies. 
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